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ABSTRACT 

The collective mode spectrum of the composite fermion state at 1/3 filling factor is 
evaluated. At zero momentum, the result coincides with the cyclotron energy at the 
external magnetic field value, and not at the effective field, in spite of the fact that only 
the former enter in the equations, thus, the Kohn theorem is satisfied. Unexpectedly, in 
place of a magneto-roton minimum, the collective mode gets a threshold indicating the 
instability of the mean field composite fermion state under the formation of crystalline 
structures. However, the question about if this outcome only appears within the mean 
field approximation should be further considered. 
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1 Introduction 



In previous papers we have introduced a generating functional approach to the 

composite fermion model of FQHE of Jain ( ||, |1) through considering the statistical 
interaction with a slowly varying parameter taken in the fermion limit 0]. In [|IJ the 
Dyson equation were solved in the Hartree-Fock approximation. The paper [0] was de- 
voted to construct a perturbative approach based on the generating functionals and mean 
statistical field values. 

In this work we continue the analysis of the implications of the general approach in 
0. Concretely the Bethe-Salpeter equation and its contribution to the linear response of 
the FQHE regimen is considered in order to discuss the elementary excitation spectrum 
predicted by the model. Technically, the discussion in this paper is very close the one in 
the work || in which the fractional statistics gas was investigated. Here, as the FQHE 
problem is examined the discussion includes the Coulomb interaction and the presence of 
the external magnetic field. The intra Landau level collective mode spectrum was already 
investigated in following an approach proposed by Feynman. A nonvanishing energy 
gap associated to magneto-roton excitations was obtained. 

Nowadays, the general point of view is that the gap for the v = 1/3 state should 
be near 0.1 e 2 /r Q , a value which is compatible with recent experimental results [fF|. 
Therefore, taking into account the also widespread interest created by the composite 
fermion description of the FQHE, in this work it is intended to determine the properties for 
the collective mode spectrum in this approach. The method employed was the numerical 
solution of the Bethe-Salpeter equation for the four-points Green's function taken in the 
first approximation for the interaction kernel. Then, as the electromagnetic response 
kernel is linearly given in terms of the four-point function at coinciding points, it follows 
H that the collective mode spectrum can also be obtained as a byproduct. 

Section 2 is devoted to review the previous introduced general background. The scheme 
of calculation is presented in Section 3. In Section 4 the results for the collective mode 
spectrum are given. 

2 Review 

In a functional way of description for the composite fermion model was proposed. 
The treatment was based in the Green's function generating functional 

Z[J h J ,r)*,r),ji,jo] = J Va i Va 'Dil)*V'ipZ[tp,ip*,A u ,ai,ao,C) 

.exp J (—ip*rj + Tj*ip + jjOj + jo^o + <j()d 2 xdx^. 

exp J (JiAi + J A ]d 3 xdx4 : ), (1) 

where the arguments of Z are sources for all the fields to be specified below and 

Z[if>,il>*,A ll ,ai,ao,C] = exp(S), ( 2 ) 

in which the action in terms of the composite fermion, statistical and electromagnetic 
fields has the form 
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the interaction vertex Vi is defined by 



Vj = — - J [i[)*(x)i/;*(x')U{x — x.')ip(x') , ip(x)} X4=x ' d d 2 xd 2 x'dx 



(3) 
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+ [if)*(x)ijj(x)U(yi — x.')} X4=x > n d xd x'dx 4 
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d x'd xdx 4 [%l)*(x)ip(x')A (x — x')ip(x')ip(x)} :i 



2m* c 2 J L v ; v yjx4-* 4 



Above, n represents the compensating charge density of the jellium, ip* and ip are 
the composite fermion fields, a M is the statistical field, A T = A e ^ + a fM + is an auxiliary 

variable and A^ is the vector potential of the external magnetic field in the gauge (A e = 
1/2B x x, At — 0), Afj, is an external electromagnetic field. The statistical field is related 
to the fermion density through the usual formula 



diix) = I Ai(x-x')tlj*(x')tlj(x')d 2 x, 
a 4 (x) = 0, 

where, the centered at x' solenoidal vector potential A is given by 



A- fx - x ) = e lJ 



ire (| x — x' |) 2 ' 
e 12 = -e 21 = l,e u =e 22 = 0, (4) 

the fractional statistical parameter 9 defines the solenoid flux in (6), which will selected 
here as equal to two flux quanta and flowing in the direction opposed to the sense in 
which the magnetic field points. This selection defines the 1/3 filling factor state M. 

In references [0-0 the Dyson equation in the Hartree-Fock approximation was solved 
and various general properties of the exact propagator and dielectric tensor obtained. One 
interesting fact following was that the mean field solutions are exactly characterized by a 
constant effective field which correspond to the external magnetic field plus an addition 
created by the statistical field. The selfenergy spectrum was evaluated and will be used 
here as a necessary piece in discussing the BS equation. 
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3 Scheme of Calculation 



In this section the technique for the evaluation of the collective mode spectrum is con- 
sidered. The general formulations developed in || were closely followed. From them, 
it becomes clear that the collective mode spectrum should appear as included in the 
spectrum of the four-point Green function. This occurs because, at least in the approxi- 
mation to be considered, the electromagnetic linear response kernel is also linearly related 
with the mentioned Green function through magnitudes which are not singular in the 
frequency. The approach followed consisted in reducing the Bethe-Salpeter equation to a 
matrix equation by using the magneto-exciton two-particle wavefunctions || and then 
solve for the frequency numerically by selecting a finite number of functions in the basis. 

The integral Bethe-Salpeter equation for the four-point Green function, after intro- 
ducing the compact notation of representing spacetime arguments by integer numbers and 
space vectors by the same numbers with arrows, can be written as 

^(1,1, 2, 2') = F {1,1', 2, 2') (5) 
J J JdAd5 d6.F (l,l', 3, 4)1^(3,4,5,6)^(5, 6, 2', 2), 



+ 



where JF Q is the interaction free four-point function and the interaction kernel in the first 
approximation is taken as the following functional derivative 

5G{5,6) 

in which W HF (1,2) is the Hartree-Fock inverse propagator 0, ||. 

As mentioned above, it is possible to simplify the integral equation (7) introducing 
magneto-exciton wave functions ||, §| which have the form 



na 



:-D n i 



d 1 
2^ z 1 

dz\ 2 



2 dz 2 



L (27r2 n + n 'n\n'\y/2 
exp[-(\z 1 \ 2 + \z 2 \ 2 + \z a \ 2 )/4}.exp[{z* 1 z 2 + z* 1 z a -z 2 zl)/2}, (6) 
z = x + yi. (7) 

These two-particle states are constructed from the one-body wave functions for elec- 
trons and holes in the mth and nth Landau levels. They are characterized by a center 
of mass momentum which can be written in complex form as z a = iq x — q y . The main 
simplification introduced by these states is that the z a is conserved and all the matrix 
elements are diagonal in this quantum number. 

Then, the integral equation (7) can be written as a linear matrix equation in the form 

•F'i nm,n' ,m' > ^) ^nn'^mm 1 ^ onm (8) 

+ J2^onm{uj) < nma\W\kla > T k i,n'm>{pi,u). 
kl 

The matrix elements of the free four point function T and the interaction kernel W 
are given by the expressions 
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( ™' 1^1 Jj' ) =JJJJ *Z'(l,4)^ (l,4|2,3)*l(2,3)dld2d3dl, (9) 

( ™' |W| \ =//// ^(l,4)W(i;4|2,3)<(2,3)dlrf2d3d4, (10) 

in which the dependence on the a quantum number has been omitted and a vertical 
representation of the bracket indices is also used. 

The interaction kernel becomes the sum of twenty contributions arising from the any- 
onic like interactions and two from the Coulomb one in the form 

20 2 



i=l 



The explicit formulae for these terms and the one for T Q are given in Appendix A. 
obtained The Feynman diagram representation of these contributions can be found in ||. 

Let us discuss now the parameters characterizing our specific problem. The statistical 
parameter will be selected as satisfying 9 = —Ikis. It corresponds with the composite 
fermions description in which an even number of flux quanta are attached to each elec- 
tron [0-0. In addition, due to the fact that we are interested in analyzing the 1/3 filling 
factor case, only one Landau will be filled, then k = 1 will be fixed M. This selection is 
related with a negative mean statistical field which cancels 2/3 of the external magnetic 
field U . All the lengths and energies here can be expressed in units of the magnetic length 
in the effective magnetic field remaining after the external field is partially compensated 
l and the cyclotron energy huj e J^ in the same field, respectively. 

The collective modes of the system, correspond with the lowest frequencies leading to 
zero eigenvalues of the inverse kernel 

J 1 " 1 = {Fo)~ l - W, (11) 

which in addition are also singularities of the dielectric response tensor component with 
both indices being temporal []5| . 

Due to the fact that T is qualitatively different for particle-hole and hole- particle 
channels it becomes useful to divide the matricial equation (10) in four sub-blocks. Con- 
cretely, the block representation is obtained by restricting, as conceived in |5j , the indices 
m and n to empty Landau levels and ml and n' to filled ones and define the block matrices 
£ and d through 

/ m! ,,,„ n' \ / m! n' \ 



and 



|W| = \€\ 

m n \ m n 



m ' \w\ n , ) = ( m ' \d\ n ' 

m n \ ... 1 



in which the properties of the magneto-exciton wave function also implies 



m , |W| n , 

m n 
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( m , \W\ n '\ = ( m ' \ d \ n ' V 
\ m n \ m n 

The explicit expressions for the matrix elements of the matrices £ and d are given in 
Appendix B. The twenty anyonic terms associated to each of these blocks S and d were 
received just by using the almost identical expressions considered in [[J. 

Now, introducing the matrix Ae as 

|Ae| ) = 5 mn 6 m r n >(e n — e n /) (12) 
m n 

in order to compact the blocks coming from the matrix representation of JF , the homo- 
geneous Bethe-Salpeter equation (10) can be reduced to the matrix form 0. 



Tiuj — Ae — £ + ir] —d 

-d + -ftu - Ae - 8* + if] 

This representation was then used to determine numerically the collective mode dis- 
persion relation. For this purpose a finite number of basis functions corresponding to the 
Landau levels index up to a maximum value were retained in constructing the BS matrix 
(15). The results are described in next section. 



(13) 



4 Results and discussion 

The roots of the determinant of the matrix (15) were found numerically after consider- 
ing only the lowest Landau level indices for the magneto-exciton basis function up to a 
maximum value n c . This number was also varied in order to check the convergence of the 
results. 

The lowest energy collective mode dispersion relation is shown in Fig.l as a function 
of center of mass momentum taken along the Xi axes. As noticed before, the energies 
are given in units of the cyclotron frequency of the effective magnetic field u e J* and the 
momentum in units of ql D . The curve corresponds to a number n c = 15 of Landau levels. 

It can be observed that the collective mode energy decreases as the center of mass 
momenta grows up to a threshold value near qlo = 2.165. At this point unstable modes 
having imaginary part of the frequency appear. These kind of modes remain up to near 
the qlo = 3.3 momentum value. Then, the unstability disappear and the real energy begin 
to grow with the momentum of center of mass. It can be noticed that unstability region 
begins for wavevector values greater than qlo — 1. 

The above described behaviour is not significantly affected by the presence of the 
Coulomb interaction. The growing and spatially periodic oscillations have a wavelength 
of the order of the magnetic length in the effective magnetic field. 

An important property which follows from Fig. 1 is that the Kohn theorem is exactly 
satisfied at zero momentum. That is, the collective mode energy became equal to three 
times Tvjj e Jl , or what is the same, equal to the cyclotron energy in the external magnetic 
field hu c . This outcome is an independing checking of the calculation in the low momen- 
tum limit. Another verification performed consisted in taking the parameter 9 = 1 by 
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also excluding the Coulomb interaction. After this, it was possible to reproduce the col- 
lective mode dispersion for a bose fluid obtained in Ref. ||. Finally, previous results for 
the inter Landau level collective excitation of Kallin and Halperin were also repeated by 
considering the coulomb interaction and assuming zero flux for the statistical solenoidal 
interaction, that is taking 9 = k = 0. 

It is known that the perturbative calculation of the energy within the composite 
fermion approach has difficulties. One of them, for example, is that after disconnect- 
ing the coulomb interaction the calculated values can hardly interpreted as reproducing 
the properties of the noninteracting electron system. Therefore, it cannot be discarded 
that the detected unstability could disappear after including correlations in the same way 
that those correlations are expected to repair the above mentioned troubles with energy. 
However, under the acceptance that the results remain valid in the exact theory, the un- 
stability could be imagined to be related with the largely discussed question about the 
role of crystal fluctuations in the FQHE ground state. Further analysis related with these 
possibilities will be considered elsewhere. 
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Appendix A 

a. The selfenergies and their statistical and coulomb contributions for each 
Landau level of index n 

£ n = e H {n) +e c (n) 

^,) = ^ ( .Ht HggjSil ) 
^»)=» + »-t +j> (-5 + *)-t( 1 / 2+ 2^tmo-£s)- 

where n is the index of the Landau level in the effective field and the constant u c is given 
by 

Uc = (-T- / 

el a mc 

b. The matrix elements of T Q 

{Tiuj — (e n — e n /) + 277) _1 n empty, n' occupied 
{—Tiuj — (e„/ — e n ) + w/) -1 n' empty, n occupied 
otherwise. 

c. The twenty two contributions to the matrix elements of interaction kernel 

W 

In the formulae below, the ^ and ^ # represents any two states of the magnetoexciton 
basis, A is the vector potential for the effective magnetic field and A12 means the solenoid 
vector (6) with spatial arguments 1 and 2. As in || dimensionless quantities e, m, c and 
H are set equal to 1, which produces energies and lengths in units of l Q and hw e Jf 

Wf = 6 1hnJ dl J d2* B (2,2)A 12 (P 1 + 1^(3,1), 
W£ = e\imJ dl J d2y* A (2,2)A 12 (P 1 +A 1 )y B (l,3), 
W£ = -ejdl J d2V B (2, 1){A X2 {P X + Ai)V* A (2, 1), 
Wf = -6 jdlj d2V A (l,2)[A 12 (P 1 + Ai)n(U), 



\ m n M 
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Wi = 9 2 J dl j d2|A 12 | 2 ^(2,2)* B (l,l), 
W£ = -6 2 j dl j d2|A 12 | 2 ^(l,2)* B (l,2), 
W 7 A = -9 2 J dl I ' d2 j d3Ai 2 Ai 3 n(l,l)^(2,3)* B (2,3), 
W 8 A = -9 2 Jdljd2j d3A 12 A 13 U(l, 3)^(1, 3)* B (2, 2), 
W 9 A = -6 2 Jdljd2j d3A 12 A 13 Tl(3, 1)^(2, 2)* B (1, 3), 
>< = -Q 2 j dl J d2 J d3Ai 2 Ai 3 n(l,2)^(3,3)* B (2,l), 
W£ = -# 2 | dljd2j d3Ai 2 Ai 3 n(2,l)^(2,l)* B (3,3), 
= -6 2 J dljd2j d3A 12 A 13 n(3,2)^(l,l)* B (2,3), 
= -9 2 j dljd2j d3A 12 A 13 U(2, 3)^(2,3)^(1,1), 
W& = 6 2 j dljd2j d3A 12 A 13 n(l,l)^(3,3)* B (2,2), 
W A = 9 2 j dlj d2j d3A 12 A 13 n(l,3)^(2,3)* B (2,l), 
W A =6 2 j dlj d2j d3Ai 2 Ai 3 n(3,l)^(2,l)* B (2,3), 
W A =6 2 j dlj d2j d3Ai 2 Ai 3 n(l,2)^(l,3)* B (2,3), 
W& = 6 2 J dlj dlj d3A 12 A 13 n(2,l)^(l,3)* B (l,3), 
W A =6 2 j dlj d2j d3Ai 2 Ai 3 n(3,2)^(l,2)* B (l,3), 
W 2 A = 6 2 j dlj dlj d3A 12 A 3t2 U(2, 1)^(3, 1)* B (2,1), 
= ~jdlj ^1^^(1,2)^(1,2), 

W 2 C 2 = jdljd2j d 3 ^j^(2,2)* B (l,l). 

where the projection operator in the first Landau level IT is defined as 

11(1,2) = ^-Lexp[-(N 2 + \z 2 \ 2 )/A + zlz 2 /2]. 

As all the vectors appearing are spatial ones the arrow over them has been omitted 
for avoid more cumbersome writing. 
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Appendix B 



a. Definitions of various auxiliary quantities 



7 n 7 *m 

r{,n.n) ( {)'"■" 



d(m, n) = (— 1) 



( 2 m+n m ! n !)l/2' 

6^ 



m+n a 



(2 m+n m!n!) 1 /2' 



m 

co{m, n) 



(m -/)!/!' 

1 



6 = -|zJ 2 
2 1 1 



fn — e ^2 Tj"> 



k=0 U 



p=l 



lm,n lib. lb. C ? u , 

p + m)\(p + n)\ 

b. Divergent integrals 

roc ^ 

Ji = e' b lim / rfr— — ^J (qr), 



= e~ b lim / dr- -e" r /2 , 

// -II ./() ' - "- 



r ^ _|_ ^ 
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The £ block matrix components for m, n > 



22 



£(m,n) = £ r (m, n), 



r=l 



£,(712,71) = -0A_L^ e -*> 

£ 2 (m, n) = -9 — - — e I - - 1 I , 



£ 3 (m,n) = —9 e , 



£4(171, n) 



\n>m 



5 mn + c(m,n) ( — f m e b + — [1 - f n -i] 



£<j(m,n) = 6 2 c(m,n)Ii, 



2 c(m,n) 

£ e (7n,n)\ n > m = -6> — - — 



(2/ 2 + X + 5m _ 1 -/i n _ 1 ), 



£ 7 (m,n)\ n > m = 9 2 



-b n—1 



E b - E r + — 



E 

k=p+l 



c(m, n) 



(1- 



1 fmlfm-! n\ 



(n — m) \ b m b n 



£ 8 (m,n) = 9 2 ^-^e- b (j- I). 



£ 9 (m,n) = 9 2 ^-^e- b 



£ 10 (m, n) = -9 2 1 j J - e b { — + 



b n + lj' 



£n (m,n) 



e2 c(m, n) _ b 



1 1 
e " [ - + 



b m + 1/ ' 



£ 12 (m,n) = -—c(m,n) 



-b n 
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S 13 (m,n) = -—c(m,n) 



<''-*> + -rE* 

k=l 



£u(m, n) = 9' 



,c(m,n) _ b l 



e -, 



2 b J 



£ 15 (m,n) = 9 2 V ' ) e~\ 



£ 16 (m,n) = 9 2 1 ' ' e~\ 



£l7Kn)| n > m = — ^ ^-^[1 - <W - /n-l] + ^l^fm j , 



£i 8 (m,n)| n > m = -— c(m,n) 



m o m m + lo" 



# 2 / m'n' A 

£i 9 (m,n)| n > m = —c{m,n) \X + g m ^ - h n _! + hm + n+l f m { 1 ~ fn)e J , 



q2 / 



£ 21 (m,n) = -u c c(m,n)e- b co ( m ' l)co(n, k)T(l + k + l + \l - k\/2){l/2f 

/=o fc=0 



.i-fc i^Ul + fc + j + I* - k\/2, 1 + \k - fj, k a | 2 /2) 

; 1 al 2V2r(i + \i-k\) 



£22(111,, n) = u c c(m,n)- 



,-b 



The d block matrix components for m, n > 



22 



d(m, n) = d r (m, n), 



r=l 
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A ( \ a d {m,n) -b( m ,\ 
di{m, n) = -9 e \^—-lJ, 

a ( \ n d ( m i n ) -b( n A 
d 2 (m, n) = -9 e y- - 1J , 

, , v Q d{m,n) ( _ b (m + n- 1)! 
d 3 (m,n) = -0 le - m — [1 - f n+m . 

d(m,n) f _ b (m + n-l)\ 
d 4 {m,n) = -9 le -n — [1 - / n+m _ 

^(m, n) = 9 2 d(m,n)Ii, 
d 6 (m, n) = (2/ 2 + X - h m+n ^) , 

j i \ n2d{m,n) f(m + n - 1)! A 

= '^-^ ^ fem+w [1 " fm + n-l}) 

4 (m,w ) = -0 2 v ' e ( — — — + t ) , 
4 \m + 1 6/ 

, / \ n 2 d ( m ' n ) -b( n i\ 

d 9 (m,n) = 2 - 7 e 6 - lj , 
i / \ n2d(m,n) _ h (\ 1 \ 

fch.) = ^(!-i), 

d 12 (m,n) = -O 2 ^^- (^(/i - h) + 2 e E > 
d 13 (m, n) = ^ 2 ^^ (Ui - / 2 ) + 2 e E > 



du(m,n) = 6 — , 



2 d(m,n) ( 1 _ b (to + ti-1)' , A 



, t \ n2 d ( m ^ n ) ( -b (m + ra-l)' A 
d w (m,n) = 9 2 4 (e ~^ &m+w [1 - /m+n-i] 1 , 



, , x n2 d(m,n) ( 1 , (to + 71 — 1)! M „ A 
dirKn) = l? 2 -^ f + &m+w + [l - iW„-i]J , 



2 d(m,n) ( _ b (m + ra-1)' , A 
d 18 (m,n) = 2 -^-^le b -n v &m+n [1 - / m+ n-i] j , 



di 9 (m, n) = f^™ft (x - hm+n ^ + e- b £ ^ , 



d 20 (m, n) = 6*^ft (x - h m+n ^ + e~ 6 ± ±) , 



ra+n 

d 21 (m,n) = -w c d(m,n)e- 6 ]T co(to, Z)(l/4)'/ 2 |z a | , / 2 r(Z + 1/2) 
1 F 1 (/ + 1/2,1 + / + 1,1^/2) 



v^(l/2)T(Z + l) 
d22ijn,n) = u c d(m,n) 
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Figure Caption 

Fig. 1 The collective mode energy in units of hw e Jf as a function of the center of mass 
momentum. Note that the Kohn theorem is satisfied at zero momentum and that the 
system develops a threshold after which unstable modes appear in a band of momentum 
values. 
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